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We study two large classes of alternative theories, modifying the action through algebraic, 
quadratic curvature invariants coupled to scalar fields. We find one class that admits solutions 
that solve the vacuum Einstein equations and another that does not. In the latter, we find a de- 
formation to the Schwarzschild metric that solves the modified field equations in the small coupling 
approximation. We calculate the event horizon shift, the innermost stable circular orbit shift, and 
corrections to gravitational waves, mapping them to the parametrized post-Einsteinian framework. 
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I. INTRODUCTION 

Although black holes (BHs) are one of the most strik- 
ing predictions of General Relativity (GR), they remain 
one of its least tested concepts. Electromagnetic ob- 
servations have allowed us to infer their existence, but 
direct evidence of their non-linear gravitational struc- 
ture remains elusive. In the next decade, data from 
very long-baseline interferometry [l|, Q and gravitational 
wave (GW) detectors [3l-[20j should allow us to image 
and study BHs in detail. Such observations will test GR 
in the dynamical, non-linear or strong-field regime, pre- 
cisely where tests are currently lacking. 

Testing strong-field gravity features of GR is of utmost 
importance to physics and astrophysics as a whole. This 
is because the particular form of BH solutions, such as 
the Schwarzschild and Kerr metrics, enter many calcu- 
lations, including accretion disk structure, gravitational 
lensing, cosmology and GW theory. The discovery that 
these metric solutions do not accurately represent real 
BHs could indicate a strong-field departure from GR with 
deep implications to fundamental theory. 

Such tests require parametrizing deviations from 
Schwarzschild or Kerr. One such parameterization at the 
level of the metric is that of humpy BHs j2l| - |23j . while 
another at the level of the GW observable is the param- 
eterized post-Einsteinian (ppE) framework (2^ . [25| . In 
both cases, such parameterizations are greatly benefited 
from knowledge of specific non-GR solutions, but few, 
4D, analytic ones are known that represent regular BHs 
(except perhaps in dynamical Chern-Simons (CS) grav- 
ity [23 [271 and Einstein-Dilaton-Gauss-Bonnet (EDGB) 
gravity j28l - l32j ). 

Most non-GR BH solutions are known through numer- 
ical studies. In this approach, one chooses a particular al- 
ternative theory, constructs the modified field equations 
and then postulates a metric ansatz with arbitrary func- 
tions. One then derives differential equations for such 
arbitrary functions that are then solved and studied nu- 
merically. Such an approach was used, for example, to 
study BHs in EDGB gravity (H-Hl. 

Another approach is to find non-GR BH solutions 
analytically through approximation methods. In this 
scheme, one follows the same route as in the numeri- 
cal approach, except that the differential equations for 



the arbitrary functions are solved analytically through 
the aid of approximation methods, for example by ex- 
panding in (a dimensionless function of) the coupling 
constants of the theory. Such a small- coupling approx- 
imation (26l . Issl . [3^ treats the alternative theory as an 
effective and approximate model that allows for small GR 
deformations. This approach has been used to find an 
analytic, slowly-rotating BH solution in dynamical CS 
modified gravity [H ■ 

But not all BH solutions outside of GR must nec- 
essarily be different from standard GR ones. In fact, 
there exists many modified gravity theories where the 
Kerr metric remains a solution. This was the topic stud- 
ied in [ssj . where it was explicitly shown that the Kerr 
metric is also a solution of certain f{R) theories, non- 
dynamical quadratic gravity theories, and certain vector- 
tensor gravity theories. Based on these fairly generic ex- 
amples, it was then inferred that the astrophysical obser- 
vational verification of the Kerr metric could not distin- 
guish between GR and alternative theories of gravity. 

Such an inference, however, is not valid, as it was later 
explicitly shown in (26| . Indeed, there are alternative 
gravity theories, such as dynamical CS modified gravity, 
where the Kerr metric is not a solution. This prompted 
us to study what class of modified gravity theories admit 
Kerr and which do not. We begin by considering the 
most general quadratic gravity theory with dynamical 
couplings, as this is strongly motivated by low-energy 
effective string actions [Sa - liOj . When the couplings are 
static, we recover the results of (ssj . while when they are 
dynamic we find that the Kerr metric is not a solution. 
In the latter case, we find how the Schwarzschild metric 
must be modified to satisfy the corrected field equations. 
We explicitly compute the shift in the location of the 
event horizon and innermost stable circular orbit. 

Such modifications to the BH nature of the spacetime 
induce corrections to the waveforms generated by binary 
inspirals. We compute such modifications and show that 
they are of so-called second post-Newtonian (PN) order, 
i.e. they correct the GR result at 0{v'^) relative to the 
leading-order Newtonian term, where v is the orbital ve- 
locity. We further show that one can map such correc- 
tions to the parameterized post-Einsteinian (ppE) frame- 
work jl^ . which proposes a model- independent, wave- 
form family that interpolates between GR and non-GR 
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waveform predictions. This result supports the sugges- 
tion that the ppE scheme can handle a large class of 
modified gravity models. 

The remainder of this paper is organized as follows. 
Sec In] defines the set of theories we will investigate and 
computes the modified field equations. Sec. IIIII solves 
for BH solutions in this class of theories. Sec. IIVI dis- 
cusses properties of the solution and Sec. |V] studies the 
impact that such BH modifications will have on the GW 
observable. Sec. IVII concludes by pointing to future pos- 
sible research directions. For the remainder of this pa- 
per, we use the following conventions: latin letters in 
index lists stand for spacetime indices; parentheses and 
brackets in index lists stand for symmetrization and an- 
tisymmetrization respectively, i.e. A(^ab) = (^a6 + ^&a)/2 
and A 
G 



ab] = {Aab — ^&a)/2; wc usc gcometric units with 
c = 1. 



II. QUADRATIC GRAVITY 

Consider the wide class of alternative theories of grav- 
ity in 4-dimensions defined by modifying the Einstein- 
Hilbert action through all possible quadratic, algebraic 
curvature scalars, multiplied by constant or non-constant 
couplings: 

S = J d^x^ {kR + aih{d)R' + a2f2{^)RabR'''' 



a^h{^)RabcdR''''"' + aiUi^) Rated* R" 
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where g is the determinant of the metric gat, 
(R, Rab, Rabcd,* Rabcd) arc thc Ricci scalar and tensor, 
the Riemann tensor and its dual 113 respectively, £inat 
is the Lagrangian density for other matter, ■)? is a scalar 
field, {ai,(3) are coupling constants and k = (IGttG)"^. 
All other quadratic curvature terms are linearly depen- 
dent, e.g. the Weyl tensor squared. Theories of this type 
are motivated from fundamental physics, such as in low- 
energy expansions of string theory [37H40| ■ 

Let us distinguish between two different types of the- 
ories: non-dynamical and dynamical. In the former, all 
the couplings are constant (//(i?) = 0) and there is no 
scalar field {(3 — 0). Varying Eq. ([T]) with respect to thc 
metric and setting fi{'d) = 1, we find the modified field 
equations 
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where T"^"' is the stress-energy of matter and 

■Hab = 2RabR - \gabR^ - 2V abR + 2gabOR , (3a) 
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with Vq, Vab ~ VaVb and □ = VaV" thc first and 
second covariant derivatives and the d'Alcmbcrtian, and 
using the Weyl identity ACa^'^Cbcde = gabCcdefC^''^^ , 
with Cabcd the Weyl tensor. 

The dynamical theory is specified through the action 
in Eq. (jT]) with fi{^) some function of the dynamical 
scalar field with potential V{i!)). For simplicity, we re- 
strict attention here to functions that admit the Taylor 
expansion fi{^) = /,(0) + /;(0)i? + 0{^'^) about small i?, 
where /i(0) and //(O) are constants. The i9-indcpendcnt 
terms, proportional to /i(0), lead to the non-dynamical 
theory, and we thus ignore them henceforth. Let us then 
concentrate on /i(??) = qi?, where we reabsorb the con- 
stants Ci = f[ into oii, such that Oiifi{'d) — > aid. The 
field equations are then 
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where T^;> = f [V^i^V^i? 
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the scalar field stress-energy tensor and 
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- 2i?,(,VSb) + gab (I'^Vei? + R'^'VcVd) 
+ d [2R"'Racbd - ^abR + aRab 
+ ^gabiOR-RcdR''')] , 

VcRab) + ARacbdV'v"^ 
1? [2 {RabR ~ AR'^'Racbd + Vahi? - 2URab) 

\gab {R^ - ARcdR'% , 
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with Va = Vai? and the Levi-Civita tensor. Notice 
that a^lCab = OLcsGab, where acs and Gab are the CS cou- 
pling constant and the CS C-tensor [231 ■ The dynamical 
quadratic theory includes dynamical CS gravity as a spe- 
cial case. Variation of the action with respect to i? yields 
the scalar field equation of motion 

pUd - = -aiR^ - a2RabR'''' 



CtsRabcdR 



abed 



0:4Rabcd*R 



abed 



(6) 



Both thc non-dynamical and dynamical theories arise 
from a diffeomorphism invariant action, and thus, they 
lead to field equations that are covariantly conserved, 
i.e. the covariant divergence of Eq. ([2]) identically van- 
ishes, while that of Eq. ([5]) vanishes upon irnpqsition of 
Eq. ([5]), unlike in non-dynamical CS gravity |27| . 
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III. NON-SPINNING BLACK HOLE SOLUTION 
A. Non-dynamical Theories 

The modified field equations of the non-dynamical the- 
ory have the interesting property that metrics for which 
the Ricci tensor vanishes are automatically solutions. 
One can see that if Rab = 0, then Eqs. (|3a |) - (pc| van- 
ish exactly, thus satisfying the modified field equations 
in Eq. This generalizes the result in |35|, as we here 
considered a more general action. 

The reason for this simplification is the Gauss-Bonnet 
and Pontryagin identities. The integral of the Gauss- 
Bonnet term G = ~ ^RabR"'' + RabcdR"'""^ is propor- 
tional to the Euler characteristic £, while that of the Pon- 
tryagin density Rabcd* R'^'^'^'^ is proportional to the Chern 
number C. Thus, the RabcdR"'""^ and the Rabcd* R"'""^ 
terms can be removed from the action in Eq. ([l} in favor 
of S and C. Since the variation of these constants van- 
ishes identically, the field equations can be rewritten to 
depend only on the Ricci tensor and its trace. 

This feature has a natural generalization for a wider 
class of alternative theories of gravity. If an action for 
an alternative theory contains the Riemann tensor or 
its dual only in a form that can be rewritten in terms 
of topological invariants (with no dynamical couplings), 
then the field equations will be free of Riemann, and 
thus, all vacuum GR solutions will also be solutions of 
such modified theories. Therefore, any action built from 
powers of the Ricci scalar or products of the Ricci tensor, 
possibly coupled to dynamical fields, and with Riemann 
tensors entering only as above, admits all vacuum GR 
solutions. 

These results have important consequences for at- 
tempts to test GR in the strong field. Electromagnetic 
GR tests that aim at probing the Kerr nature of BHs 
would be insensitive to such modified theories. On the 
other hand, observations that probe the dynamics of the 
background, such as GW observations would be 

able to constrain them. 



B. Dynamical Theories 

The modified field equations in the dynamical the- 
ory, however, are not as simple, as clearly they are 
not satisfied when Rab = 0. This is because ^/^^ de- 
pends on V^v'^Rabcd and /C^^^ depends on V^Wc *R{a^b)'^- 
Let us search for small deformations away from the 
GR Schwarzschild metric that preserve stationarity and 
spherical symmetry. The only relevant term here then is 
"^ib ' ^ ^ab vanishes in spherical symmetry, as already 
analyzed in [2^. 

We thus pose the ansatz 



-fo [1 + ehoir)] dt^ + fo' [1 + efco(r)] dr^ 



and -d = d + e-3, where /o = 1 - 2Mo/r, with Mq the 
"bare" or GR BH mass and (t, r, 9, 0) are Schwarzschild 
coordinates, while dU^ is the line element on the 2-sphere. 
The free functions (/io,fco) are small deformations from 
the Schwarzschild metric, controlled by a function of the 
coupling constants {ai,f3) that we define below; e is a 
book-keeping parameter. 

Before we solve the field equations, let us discuss 
the scalar field potential V{i}). There are two distinct 
choices for this potential: a flat (V'{'&) = 0) or non-flat 
{V'{'d) 0) potential. For the non-flat case, the poten- 
tial must be bounded from below for the theory to be 
globally stable, and thus it will contain one or more min- 
ima. The scalar field would tend towards the minimum 
of the potential, where the latter could be expanded as a 
quadratic function about the minimum (assumed here to 
be at zero): V « ^m'^'d'^. One might treat the flat poten- 
tial as the limit — >■ of the above non-fiat potential, 
but this limit is not continuous at the point = 0. 
The massive case must thus be treated generically and it 
turns out to be sufficiently complicated that we restrict 
our attention only to the massless (flat) casef4l|. 

With this ansatz, we can solve the modifled field equa- 
tions and the scalar field's equation of motion order by 
order in e. Through the small-coupling approximation, 
we treat a — 0{e) and /3 = 0{e). To zcroth-order in 
e, the field equations are automatically satisfied because 
the Schwarzschild metric has vanishing Ricci tensor. To 
this order, the scalar field equation can be solved to find 
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This is the same solution found in [36| for dilaton hair 
sourced in EDGB gravity. The scalar field depends only 
on as , since the term proportional to 014 vanishes identi- 
cally in a spherically symmetric background. 

We can use this scalar field solution to solve the mod- 
ified field equations to 0{e). Requiring that the metric 
be asymptotically flat and regular at r = 2Mo, we find 
the unique solution ho = T{l + ho) and fco = —T{l + ko), 
where J" = -(49/40) C (Mo/r) and 
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and where we have defined the dimensionless coupling 
function ( = al/{l3KM^) ~ 0[t). This solution is the 
same as that found in EDGB gravity [42|. Our anal- 
ysis shows that such a solution is the most general for 
all dynamical, algebraic, quadratic gravity theories, in 
spherical symmetry. 

The demand that the metric deformation be regu- 
lar everywhere outside the horizon has led to a term 
that changes the Schwarzschild BH mass, i.e. there is 
a correction to gu and g„ that decays as l/r at spa- 
tial infinity. We can then define the physical mass M = 
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Mq [1 + (49/80) C] , such that the only modified metric com- 
ponents become gu = —/(I + h) and g^r = /~^(1 + k) 
where h = C/(3/)(A//r)3/i and k -{(,/ f){M /rfk, and 
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and where / = 1 — 2M/r. Physical obscrvablcs arc re- 
lated on the renormalized mass, not the bare mass. This 
renormalization was not performed by (42j . 

In fact, one need not fix the single constant of integra- 
tion which appears in finding this solution. Any value 
of the integration constant, after renormalization, is ab- 
sorbed into the renormalized mass. Rather than a family 
of spacetimes, there is a unique spacetime after renor- 
malization. 

The sign of the coupling constant can be determined 
by computing the energy carried by the scalar field in 
Eq. dH]). The energy is E(^^) = j^T^J^hH^-i^/'^cPx, where 
E is a i = const, hypersurface outside of the horizon (so 
that it is spacelike everywhere), = [d/dt]"^ and 7 is the 
determinant of the metric intrinsic to E. Wc find that 
E[^) = {9/7)(^kttM. For stability reasons, we require 
that > 0, which then implies C ^ and a|//3 > 0. 

Although we here considered non-spinning BHs, our 
analysis can be generalized to spinning ones, by sepa- 
rating the theory and its solutions into parity-even and 
parity-odd sectors. A parity transformation consists of 
the reficction —x^, which for a spinning BH metric 
implies a — > —a, where \S^\ = M\a\ is the magnitude 
of the spin angular momentum. Expanding the spinning 
BH solution as a power series in a/M, we see that the 
Kretschmann scalar Rabcd.R"'^'^'^ has only even powers of 
a/M (even parity sector), while the Pontryagin density 
*RR has only odd powers of a/M (odd parity sector). 
These quantities source the "d equation of motion, there- 
fore driving even and odd metric perturbations respec- 
tively. The solution found here is of even parity and cor- 
responds to the 0(0'') part of the metric expansion for a 
slowly-spinning BH in dynamical quadratic gravity. The 
next order, O(a^), is parity odd and is sourced only by 
the Pontryagin density, since i?^, RabR"^^ , and Rabcd.R"'^'^'^ 
arc all even under parity. The solution sourced by just 
the Pontryagin density is identical to that in dynamical 
Chern-Simons gravity (all 0;^ = except for 04) and was 
found in [2^ . From the parity arguments presented here, 
we see that the exact same modification arises at 0{a^) in 
the more general dynamical quadratic gravity considered 
here. Therefore, to O(a^), the modification in dynamical 
quadratic gravity is simply the linear combination of the 
OjaP) solution found here and the 0{a^) solution found 
in [2I. 



IV. PROPERTIES OF THE SOLUTION 

The solution found is spherically symmetric, station- 
ary, asymptotically flat, and regular everywhere except 
at r = 0. It represents a non-spinning BH with a 
real singularity at the origin, as evidenced by calculat- 
ing the Kretschmann scalar expanded to 0(0 : K = 
RabcdR'''""^ = K ~ 32CM^/r'^K, where K = A8M'^/r^ 
and 
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The location of the event horizon, i.e. the surface of infi- 
nite redshift, can be computed by solving gu = to find 
Veh/M = 2 — (49/40)C- The metric remains Lorentzian 
(i.e. sgii{g) < 0) everywhere outside r^n provided C is 
sufficiently small (specifically, < ( < (120/361)). 

One can also study point-particle motion in this back- 
ground. Neglecting internal structure and spins, test- 
particle motion remains geodesic (isj and the equation 
of motion reduces to 7-^/2 = V^^"^ -I- SV^a, where the over- 
head dot stands for differentiation with respect to proper 
time and 
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(13) 

where (E, L) are the conserved quantities induced by the 
timelikc and azimuthal Killing vectors, i.e. the particle's 
energy and angular momentum per unit mass. 

One can solve for the energy and angular momentum 
for circular orbits (43 | through the conditions r = and 
l/'f = to find E = i?GR -I- 5E and L = L^r + 6L, where 
Eou = /(I - 3A//r)-i/2, = {Mrf'^E^^Jf and 
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From this expression, we can find the modified Kepler 
law by expanding a; = L/r^ in the far field limit: 
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where u)"^^ 



M/r^[l + 0{M/r)]. If in addition to 
the above circular orbit conditions one evaluates the 
marginal stability condition 1/,". = 0, one finds that the 
shift in the ISCO location is 
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V. IMPACT ON BINARY INSPIRAL GWS 

As evidenced above, such a modified theory will intro- 
duce corrections to the binding energy of binary systems. 
Consider a binary with component masses mi 2 and to- 
tal mass m ~ mi + m2. The binding energy, to leading 
0{m/r, 0, can be obtained from iScn and SE in Eq. 
by the transformation mim2 — >■ rri^r] and expanding in 
M/r ^ 1. This trick works to leading order in C, and in 
m/r only and it leads to 



^ + 1(7 
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Using the modified Kcplerian relation of the previous Sec- 
tion, this becomes 

Eb{F) = - i {2TrmFf - ^mr^C i2nmFf , (19) 

to leading 0{mF,(^)^ where F is the orbital frequency 
and 77 = mim2/m? is the symmetric mass ratio. Such a 
modification to the binding energy will introduce correc- 
tions to the binary's orbital phase evolution at leading, 
Newtonian order. 

A calculation of the phase and amplitude waveform 
correction that accounts only for the leading-order bind- 
ing energy modification is incomplete. First, higher 
0{m/r) terms in Ei, are necessary for detailed GW tests. 
These terms, however, are not necessary to find the 
leading-order, functional form of the waveform correc- 
tion; this is all one needs to map these modifications to 
the ppE scheme. 

To be consistent, we must also consider the energy fiux 
carried by the scalar field. This program involves solv- 
ing for the perturbation on top of the background solu- 
tion given in Eq. dH). The solution can be found using 
post-Newtonian integration techniques and is in prepa- 
ration [45| . The modification to radiation reaction due 
to the scalar field is subdominant (of much higher post- 
Newtonian order) compared to the modification to the 
binding energy calculated here, as will be shown in a 
forthcoming paper (45j . 

Let us now compute the orbital phase correction due 
to modifications to the binding energy. The orbital phase 
for a binary in a circular orbit is simply 



0(F) 



{E') {E)-^ujduj, 



(20) 



where uj = 2'kF is the orbital angular frequency, E' = 
dE/dw and E = —{'i2/b)rfm?r'^Lo^ is the loss of binding 
energy due to radiation. This expression for E is the GR 
quadrupole form, which was shown [l7| to be valid in 
the small-coupling limit in asymptotically flat spacetimes 
when the action is of the form we use. Neglecting E^'^'^ 
and to leading 0{muj, (), the orbital phase 
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where the GR phase is (/)gr = — 1/(327]) (27rmF)~^/^. 
The leading order correction is of so-called 2PN order, 
as it scales with {m F)'^/^ (down by relative to the 

leading-order GR result. 

Similarly, we can compute the correction to the 
frequency-domain GW phase in the SPA, by assuming 
that its rate of change is much more rapid than the GW 
amplitude's. This phase is (see e.g. [46| ') 



20(to) - 27r/to , 



(22) 



where to satisfies the stationary phase condition F{ta) = 
f 12. with / the GW frequency. Neglecting E^^^ and to 
leading 0{mLd, Q, we find that 
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where u = irAif is the reduced frequency and M. = 
jy3/5j^ is the chirp mass. Similarly, the Fourier- domain 
amplitude scales as \h\ cx F(to)~^^^, which then leads to 
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where |/igr| is the GW amplitude in GR. In principle, 
there could be additional corrections to \h\ from modi- 
fications to the first order equations of motion, but 
has shown that these vanish in the small coupling ap- 
proximation. 

The modifications introduced to the inspiral waveforms 
can be mapped to parametrized waveform models that fa- 
cilitate GR tests. In the ppE framework , the simplest 
parameterization is 

h = \h\aAl + av'W^) cxp[i*^^(l + Pv^n")] , (25) 

where (a, a, /3, 6, c, d) are ppE parameters. Our results 
clearly map to this parameterization with a = (5/6)C, 
/3 = (50/3)C, a = 4/3 = 6 and c = -4/5 = d. Since 
the radiation carried by the scalar field is of higher post- 
Newtonian order, including it will not change these ppE 
parameters. Future GW constraint on these parameters 
could be translated into a bound on the class of alterna- 
tive theories considered here. 

Preliminary studies suggest that GW detectors, such 
as LIGO, could place interesting constraints on the pa- 
rameter /3. Given a signal-to-noise ratio of 20 for a com- 
parable mass binary inspiral signal, one might be able to 
constrain /3 < 10""'^ when 6 = 4/3 [43|. This bound would 
translate to a ^-constraint of C < 10~^, which should be 
compared to the current double binary pulsar constraint 
C ^5 10'' mi- We then see that GWs could place much 
stronger constraints on non-linear strong field deviations 
from GR relative to current binary pulsar bounds. 



VI. FUTURE WORK 

The study presented here shows that there is a wide 
class of modified gravity theories where Schwarzschild 
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and Kerr are not solutions, yet their waveform modifi- 
cations can be mapped to the ppE scheme. This study 
could be extended by investigating higher-order in v, PN 
corrections to the waveform modifications. Such a calcu- 
lation would require one to solve for the two-body metric 
in this specific class of theories. Although this can in 
principle be done within the PN scheme, in practice the 
calculation will be analytically quite difficult, due to the 
non-linear terms introduced by the modified theory. 

Another possible extension is to investigate the effect 
of different potential terms to the results presented here. 
For example, one could postulate a cosine potential and 
see how this modifies the solutions found. Such cosine 
potentials arise naturally due to non-linear interactions in 
effective string actions. The inclusion of such a potential 
will probably render the problem non-analytic, forcing 
us to solve the equations of motion for the scalar field 
numerically. 

One other avenue of future research is to find analytic, 
closed form solutions for BHs rotating arbitrarily fast 
in dynamical quadratic gravity. The analysis presented 
here applies only to non-rotating BHs, and we have dis- 
cussed how it would be modified when considering slowly 



rotating BHs. Exact, closed form solutions for rapidly 
rotating BHs, however, remain elusive. One might have 
to integrate the equations numerically to find such solu- 
tions. One possible line of attack is to evolve the field 
equations in a 3 -|- 1 decomposition, starting with a dense 
and rotating scalar field configuration. Upon evolution, 
this scalar field will collapse into a rapidly rotating BH, 
yielding a numerical representation of the solution one 
seeks. 
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